Irreversible Aggregation and Network Renormalization 



Seung-Woo Son, 1 Golnoosh Bizhani, 1 Claire Christensen, 1 Peter Grassberger, 1,2 and Maya Paczuski 1 

Complexity Science Group, University of Calgary, Calgary T2N 1N4, Canada 
2 FZ Jiilich, D-52425 Jiilich, Germany 
(Dated: March 7, 2011) 

Irreversible aggregation is revisited in view of recent work on renormalization of complex networks. 
Its scaling laws and phase transitions are related to percolation transitions seen in the latter. We 
illustrate our points by giving the complete solution for the probability to find any given state in 
an aggregation process (k + 1)X — > X, given a fixed number of unit mass particles in the initial 
state. Exactly the same probability distributions and scaling are found in one dimensional systems 
(a trivial network) and well-mixed solutions. This reveals that scaling laws found in renormalization 
of complex networks do not prove that they are self-similar. 

PACS numbers: 89.75.Hc, 02.10.Ox, 05.70.Ln 



Droplets beget rain, goblets coagulate to make butter 
or cream, and dust particles stick together to form ag- 
gregates that can eventually coalesce into planets. At 
the microscopic level, irreversible aggregation of atoms 
and molecules creates many familiar forms of matter 
such as aerosols, colloids, gels, suspensions, clusters and 
solids [1]. Almost a century ago, Smoluchowski proposed 
a theory based on rate equations to describe processes 
governed by diffusion, collision and irreversible merging 
of aggregates [2{. The theory predicts how many small 
and large clusters exist at any given time and yields a 
mass distribution that depends on certain details such 
as the initial conditions, reactions present, relative rates, 
the presence or absence of spatial structure, etc. A key 
interest to physicists has been to derive scaling laws that 
characterize different universality classes [3, and refer- 
ences therein]. 

By contrast, wide interest in complex networks 0- 
0) has emerged recently. Vast applications to physics, 
computer science, biology, and sociology and ref- 

erences therein] continue to be vigorously investigated. 
An important question is whether or not complex net- 
works exhibit self-similarity at different length scales and 
if they can be grouped into universality classes on that 
basis. Renormalization schemes for networks were pro- 
posed [n|4l4j to address this question. Scaling of the 
mass or degree distribution of the renormalized nodes 
was used to argue that many complex networks are self- 
similar. The semi-sequential renormalization group (RG) 
flow underl ying the box covering of [lll - ll4j was studied 



mg 

carefully in 15|, [16( , where it was found that scaling laws 



may be related to an "RG fixed point" which was ob- 
served for a wide variety of networks. A convenient, fully 
sequential scheme called random sequential renormaliza- 
tion (RSR) was introduced [17(. At each RSR step, one 
node is selected at random, and all nodes within a fixed 
distance £ of it are replaced by a single super-node. 

We point out a simple mapping between RSR and irre- 
versible aggregation on any graph. Hence any conclusion 
drawn for one process holds also for the other. Indeed, 



a local coarse-graining step to produce a new super-node 
represents one aggregation event, where a 'molecule' ag- 
gregates with all its neighbors within distance I to pro- 
duce a new cluster. Exact analysis in one dimension 
reveals that even this trivial network exhibits scaling 
laws for the cluster mass distribution under RSR - with 
exponents that depend on £. Consequently, and some- 
what counter-intuitively, self-similarity observed in RSR 
and similar network renormalization schemes cannot be 
used to prove that complex networks are themselves self- 
similar. Instead scaling laws arise due to a percolation 
transition in irreversible aggregation. 

The correspondence between aggregation and renor- 
malization is relevant for any model with stochastic 
coarse-graining of a network. For instance, the theory 
of space and time "Graphity" [3 EH, based on loop 
quantum gravity, involves a stochastic coarsening simi- 
lar (albeit more structured) to RSR. Hence the critical 
point of aggregation may also be relevant in that and re- 
lated cases. The breakdown of conventional universality, 
where critical exponents depend on the microscopic scale 
of coarse-graining, £, seems to present a dilemma for the- 
ories based on stochastic coarse-graining of a network to 
arrive at e.g. a universal large scale theory of gravity. 

In order to demonstrate these points, here we consider 
irreversible aggregation (k+l)X — > X, where a randomly 
picked cluster coalesces with k neighbors. For even k — 
21 this corresponds precisely to RSR on a 1-d chain with 
coarsening range I. The mass of the newly formed cluster 
is the sum of the (k + 1) masses. We assume that the 
'target' cluster is picked with uniform probab ility from 
all clusters. Other choices will be discussed in [201 ] . 

Let us start with the model defined on a ring, i.e., with 
periodic boundary conditions. Initially, Nq sites labelled 
by i 6 [1, ...N ] are each occupied by a particle of mass 
m = 1. Time can be either discrete or continuous, but 
we demand that two events never happen simultaneously. 
Hence events, ranked by increasing time, are denoted by 
positive integer values t. For each event, particles coag- 
ulate to form clusters of mass in > 1. More precisely, an 
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FIG. 1: (Color online) Illustration of aggregation on a ring 
with k = 1, No = 24, and N — 5. The tree in color corre- 
sponds to a cluster of mass m — 5. It has five leaves (blue) 
and four internal nodes (red). Its leaves start at site i and 
end at site i + m — 1. The numbers beside internal nodes 
correspond to the time when coalescence occurs. 

event consists of picking a random cluster with uniform 
probability and joining it with k clusters to its imme- 
diate right, using periodic boundary conditions. For k 
even, the same results are found if we aggregate clusters 
symmetrically. After t events, N t = No — kt clusters exist. 
Our main result is the probability to find any sequence 
of adjacent cluster masses p N °(m,i, m-i . . . mj\r t ) - where 
a cluster of mass mi is followed by a cluster of mass 7712, 
etc., moving clockwise (see Fig. [1]). we start with the 
single cluster mass probability. 

Cluster masses are restricted to m = 1 (mod k). Defin- 
ing m — 1 = ks, the integer s is the number of events 
needed to make the cluster of mass m. As depicted in 
Fig. [TJ we can represent any realization of the process 
by a forest of N t rooted trees with N leaves and t in- 
ternal nodes. Each tree a has s a internal nodes, with 
J2 a s a — t- We simplify the notation by N for N t . 

Let 7r^°(m) denote, for fixed k (the dependence on k 
is not written explicitly in the following), the probability 
that a cluster of mass m has its left-most member at site 
i G [1, No] after t events. The probability that any of the 
N clusters picked at random has mass m is then 
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because there are No choices for i and the chance to 
pick that particular cluster, given that it exists, is 1/N. 
Since events occur completely at random, each history 
occurs with equal probability. The term 'history' refers 
to a fixed forest, which includes a fixed temporal order 
of events. Thus TT^°(m) is equal to the number of histo- 
ries leading to a final configuration with a cluster of mass 
m starting at position i, divided by all possible histories 
leading to N clusters. The latter is equal to 



where each of the t factors equals the number of choices 
for the next event. Using Pochhammer k— symbols or, 
equivalently, generalized rising factorials 2lM24| , this can 
be written as rohi s t,tot = (N+k) t k . Similarly, the number 
of histories leading to a cluster of size m starting at a 
fixed position i is 

"•hist, duster = (m - k)(m - 2k) x . . . 1 = (l) s>k (3) 

and the number of histories for the remaining N — 1 clus- 
ters is 

"hist.rost = (N - m - k)(N - m - 2k) x . . . (N - 1) 
= (*-!)*_.,* ■ (4) 

So far we have not included the number of choices associ- 
ated with different time orderings for the s events in the 
cluster and (t — s) events in the rest of the forest. The 
number of different time orderings is given by 
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Combining Eqs. (JXJ) to ([S]), we obtain 

PN{ ] ~ n[s) (N + k) uk 
t\ (N-l) t _ Sik (l) s>k 
s) (AO 
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This result can be further simplified into beta functions 
or, more conveniently, fc-beta functions (see e.g. [22j), 
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giving a remarkably simple final result 

't\ Bk(N - m, m) 
B k (N -1,1) 



(7) 



We make a number of observations: (1) For k = 1 
the process maps to bond percolation on a ring. For 
N = 2, the mass distribution is uniform over the en- 
tire range m £ [l,Ao — 1]. For N > 2, the distribu- 
tion is proportional to the (N — 2) nd factorial power 
((N - m- 1)(N - m - 2) • ■ • (N - m - N + 2)). (2) 
For N = 2 and any k > 1, p^°(m) is symmetric under 
the exchange m <-> No — m. (3) For N = 2 and k = 2 
we obtain an equation formally identical to Spitzer's dis- 
crete arcsine law for fluctuations of random walks [25j |. 
(4) Asymptotic power laws for No — > oc can be deter- 
mined using Stirling's formula. If N is fixed and both m 
and (No — m) — > oo, 



n h ist,tot = A/q x (N - k) x . . . (N + k), 
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FIG. 2: (Color online) Cluster size distributions after t = 
50 events for k — 2, for different values of N averaged over 
10 6 realizations compared to exact results. The large size 
behavior changes from increasing to decreasing power law at 
N — k + 1. Inset: The discrete arcsine law found for N = 2. 



For small masses, this gives a decreasing power law with 
exponent — 1 + 1/fc. For N = k + 1, the power law 



-V fe holds up to the largest possible value, 
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m = No — N + 1, and the cutoff is a step function. For 
m/No — > 1 different power laws appear if TV ^ k + 1, 
and the sign of the exponent changes at N = k + 1. For 
TV < k + 1, the distribution has a peak at m/No — > 1, 
while it goes to zero for N > k + 1. These scaling laws 
are illustrated for k = 2 in Fig. [2] (5) The scaling laws 
found for m ^ No are identical to those obtained by 
Krapivsky [26| for the well mixed case. However, the be- 
havior for m/No — > 1 given in 126] does not agree with 
our result. (6) The probability p N ° (to) satisfies a number 
of recursion relations: 

*m + *) = , m jf^N N+1) ^ *W ■ 
(m + k — l)(iVo — to — k) 
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A third nonlinear recursion relation is given later. 

Joint distributions for masses of adjacent clusters can 
also be found. We denote by p 1 ^- (mi , to 2 ) the probability 
to find a cluster of mass mi followed immediately to the 
right by a cluster of mass m 2 . This is non-zero only if 
mi = 1 + sik and m 2 = 1 + s 2 k, where s a is the number 
of events needed to form a cluster of mass m a . By the 
same arguments that led to Eq. © we get 

V s o,si,s 2 / (N) tk 

where so = ^ — X)fl=i S P an< ^ ^ ne nrs ^ factor is the multino- 
mial coefficient instead of the binomial coefficient. When 



a = 2, it is a trinomial coefficient that counts the number 
of ways in which the three sequences of events - for the 
two clusters considered, and for all ( N — 2) other clusters 
- can be interleaved in a single history. 

For any 1 < a < N — 1 the joint probability distri- 
bution for a consecutive, adjacent clusters is a product 
of a multinomial coefficient and a + 1 Pochhammer k- 
symbols, divided by the Pochhammer /c-symbol related 
to the total number of possible histories given TVo initial 
particles. Defining again so as the number of events in all 
clusters except the first a ones, we can write the result 
compactly as 
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In particular, this can also be done for the joint distri- 
bution for all N masses by setting a = N — 1. The re- 
sulting expression is then manifestly invariant under any 
permutations of N numbers (mi, . . .mjsf). Hence the TV- 
cluster probability is independent of the spatial ordering 
of the clusters. While there are obvious correlations be- 
tween the mass values (the sum of all cluster masses must 
be iVo), there are no spatial correlations. 

We now consider a line of No particles with open 
boundaries. Again, aggregation events consist of a ran- 
dom choice of a cluster, followed by its amalgamation 
with its k nearest neighbors to the right. The target 
cluster must be at least k steps away from the right-most 
boundary. Following the same arguments leads immedi- 
ately to Eq. © for a = N — 1, showing that the two 
models lead to precisely the same statistics. 

The absence of spatial correlations indicates that the 
same dynamics might result for the well-mixed case. Now 
we start with a bucket containing No balls, each of unit 
mass. An event consists of taking k + 1 balls out of the 
bucket, merging them together, and returning the new 
ball to the bucket. The k+ 1 balls are chosen completely 
at random, independent of their masses. 

The single cluster mass distribution for the well-mixed 
model can be obtained using the same strategy as be- 
fore, but the details are quite different. Consider the 
total number of histories. Since events now correspond 
to choosing any k + 1 balls out of No — kt balls, we have, 
instead of the Pochhammer fc-symbol, a product of bino- 
mial coefficients, 
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The expressions for "hist, cluster an( i n hist,rost are analo- 
gous, with the factors (to — jk) (resp. (iVo — to — jk)) in 
Eq. ((3|) (resp. ((!])) replaced by binomial coefficients. The 
number of time orderings "ordorings is exactly the same 
as before, but the first factor No/N in Eq. (O has to 
be replaced by 4° ( N °) . Putting all these things together, 
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many cancellations take place, leading exactly to Eq. (J7J. 
This argument can be similarly extended to get the full 
TV-particle distribution function, obtaining exactly the 
same result as before, for any k. 

The time-reversed process of aggregation is fragmen- 
tation. When considering the fragmentation process as- 
sociated with any of these models, we have to care- 
fully evaluate fragmentation rates. Assuming uniform 
rates would not lead to all time-reversed histories hav- 
ing the same probability. Indeed the fraction of all 
mergers associated with making a cluster of mass m' is 
(s'/t) — (m' — I) /(No — N), which must equal the proba- 
bility that an existing cluster of mass m! will fragment at 
the next step in the time-reversed process. If it does, then 
for consistency its fragmentation products must have a 
mass distribution given by p™ +1 (m). A quadratic recur- 
sion relation for p^ +k (m) can then be obtained by con- 
sidering the likelihood of all fragmentation events in a 
configuration of N clusters, with m being the mass of 
one of the resulting k + 1 fragmentation products. The 
relation is 

No . , _ No ^J- 1 N(m> - l) P ffi(mQp^i(m) 
P N+k {m) - 2^ Nn-N 



in 
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where the prime on the summation symbol indicates that 
m! must increase in steps of k. 

In summary, we derived complete solutions for the 
probability to find any given state in three models - well- 
mixed solutions, particles on a ring reacting with their k 
nearest neighbors, and the same reaction for particles on 
a line with open boundaries - and show that these solu- 
tions are precisely the same. The fact that we could solve 
exactly a one dimensional model without detailed balance 
might seem surprising since such models are in general 
not solvable. It stems from the fact that spatial correla- 
tions, although a priori not excluded, are in fact absent. 
Related to this is our finding that the well-mixed mod- 
els have exactly the same solutions. Our method can be 
used to solve the model where the target cluster is picked 
with a probability proportional to its mass [2(j . Perhaps 
generalizations of these observations hold true for more 
complicated models, in which case weighted path inte- 
grals would replace sums over histories. 

We have pointed out a direct mapping between irre- 
versible aggregation and RSR. The latter was motivated 
by claims that one can define finite fractal dimensions for 
real networks using similar but more complicated 
and ambiguous schemes. Results for RSR with I = 1 
on various graphs (critical trees [13], Erdos-Renyi and 
Barabasi- Albert networks [27], and regular lattices 28[) 
concur with our present conclusions for k = 2. Apart 
from studying a system that is sufficiently simple to be 
exactly solvable and that is obviously not fractal, here we 
presented results for t > 1, showing that scaling laws de- 
pend in a non-trivial way on I. Results for the elementary 



network (a one dimensional line) examined analytically 
here proves that scaling under stochastic network renor- 
malization arises from an underlying percolation transi- 
tion in aggregation and does not prove fractality or self- 
similarity of the underlying graph. 

Our mapping suggests that the critical behavior of 
aggregation may also turn up in "Graphity" [HI, El 
or related models, where geometry, gravity, and matter 
emerge through an aggregation process of an underlying 
graph. "Geometrogenesis" is the complementary process 
of infinite cluster formation in irreversible aggregation. 
In that case, scaling that depends on the microscopic 
coarse-graining scale I seems to add further obstacles to 
the persistent and challenging problem to derive a large 
scale theory of gravity from microscopic graph models. 
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